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Abstract 

On the basis of the pseudo-potential plane- wave (PP-PW) method and the 
local-density- functional theory (LDFT), this paper studies energetics, stress- 
strain relation, stability and ideal strength of /3-SiC under various loading 
modes, where uniform uniaxial extension and tension, biaxial proportional ex- 
tension are considered along directions [001] and [111]- The lattice constant, 
elastic constants and moduli of equilibrium state are calculated, and the re- 
sults agree well with the experimental data. As the four Si-C bonds along 
directions [111], [Til], [111] and [111] are not the same under the loading along 
[111], internal relaxation and the corresponding internal displacements must 
be considered. We find that, at the beginning of loading, the effect of internal 
displacement through shuffle and glide plane diminishes the difference among 
the four Si-C bonds length, but will increase the difference at the subsequent 
loading, which will result in a crack nucleated on {111} shuffle plane and a 
subsequently cleavage fracture. Thus the corresponding theoretical strength 
is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. 
However, with the loading along [001], internal relaxation is not important 
for tetragonal symmetry. Elastic constants during the uniaxial tension along 
[001] are calculated. Based on the stability analysis with stiffness coefficients, 
we find that the spinodal and Born instabilities are triggered almost at the 
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same strain, which agrees with the previous molecular dynamics simulation. 
During biaxial proportional extension, stress and strength vary proportionally 
with the biaxial loading ratio at the same longitudinal strain. 
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I. INTRODUCTION 

Investigation of stability and ideal strength of materials is always an attractive issue 0" U 
due to the following facts: (1). Stability of materials is very important in elasticity theory, 
which is related with structural responses in solids, ranging from ploy-morphism, amorphiza- 
tion, and melting to fracture. 0" (2). The ideal strength of a perfect crystal represents 
an upper bound to the actual strength of crystalline materials; (3) The technology makes it 
possible to manufacture finer and finer filament, whose strength will approach the theoreti- 
cal limit. However, even for the best whisker materials. 0,0 , the realistic strength is still 
far below the predicted theoretical values. 

Recent developments of experimental technology see new opportunities of producing very 
fine nanorods (NRs) and nanotubes(NTs) & @, and the possibility of measuring their elas- 
ticity constants, strength and toughness, & §!. The typical measured diameter of NRs is 
20-30 nanometer, which can be considered free of any defects and with the ideal strength. 
Wong et al. @ measured the Young's modulus and the bending strength of silicon car- 
bon(SiC) NRs in a recent experiment. It is worth mentioning that micrometer-scale SiC 
whisker is widely used to strengthen composite materials. Thus, it is necessary to have a 
clear understanding of the stiffness, stress-strain relation, stability and strength of these 
nano sized NRs from the experimental view as well as from the theoretical point of view. 

As is well known, people are used to study the ideal strength of materials with various 
models and empirical potentials. Polanyi i! and Orowan used a model in terms of surface 
energy, interplannar space and an appropriate Young's modulus to investigate the ideal 
strength. Frenkel estimated ideal shear strength r max of a solid subjected to deformation 
of a simple shear mode. However, oversimplified functions of stress-strain were adopted in 
those methods and the functional forms for different materials were somehow arbitrary. On 
the basis of stability criteria, Milstein %' @' % investigated theoretical strength of bcc Fe, 
fee Ni and Cu with the Morse potential. The ideal strength is identified with the loss or 
exchange of stability. This method helps greatly the investigation of strength and reveals 
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a variety of interesting and surprising behaviors of materials & & H3' @; however, the 
interatomic potential they used and obtained by fitting properties of equilibrium state is 
inappropriate to be used for the investigation of the stability and strength of materials, 
which are essentially of far from being in the equilibrium state. On the other hand, the 
density functional theory(DFT) §§ , with only the input of atomic position and charge 
number , can be used to determine many structural and dynamic properties of materials 
under various conditions, including those, which are far from being in the equilibrium state. 

Through a series of comprehensive theoretical and computational studies, Hill and 
Milstein & & @ have shown that positive definiteness of internal energy is coordinate 
dependent and the stability domain depends on the choices of strain measures; while Born 
criteria are valid only under the zero load. Based on this idea, Wang et al. 0' @ 
analyzed the onset of instabilities in homogeneous lattice under critical loading and showed 
that response of the lattice is no longer a purely intrinsic property of materials, and depends 
on the applied load. Starting with these theories, Li and Wang have recently analyzed 
the stability and branching of Aluminium under various loading modes according to the first 
principle calculations. 

Heine and co-workers & §! gave a very extensive set of first principle pseudopotential 
calculations on the ploy-types in SiC. However, they gave only the bulk modulus. Lam- 
brecht et al. made a detailed investigation on the elastic constants, modulus of /3-SiC with 
the full-potential linear-muffin-tin-orbital (FP-LMTO) method. In their investigations, the 
strength of /3-SiC was obtained approximately by Orowan formula. The structural proper- 
ties of /3-SiC (ploy-type 3C) had been also investigated by the various semi-empirical models, 
e.g. , semi-empirical force model @, semi-empirical interatomic potential 0' and tight 
binding approximation & @. Modifying the Tersoff potential 0, Tang and Yip @ inves- 
tigated the lattice instability in /3-SiC and simulated the process of brittle fracture under 
hydrostatic tension based on the Hill and Milstein stability theory. The instability mode is 
the spinodal instability, and decohesion occurs as spontaneous nucleation of cracking on the 
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{111} shuffle planes. 

In the present paper, we study the energetics, the elastic constants, the stress-strain 
relation, stability and the ideal strength of /3-SiC with the density functional theory. We 
consider several loading modes, uniaxial extension and uniaxial tension along [001] and [111] 
directions, and biaxial proportional extension along [001] and [010]. The deformation is ho- 
mogeneous and elastic, and the strain can be large. The stress-strain relations are calculated, 
and the ideal strength is obtained according to the stability criteria. Owing to the unequiv- 
alence of the four Si-C bonds under the loading along [111], the internal relaxation must be 
considered and the internal displacements be calculated. With the internal displacements, 
we discuss the effect of the relaxation and failure modes. The stability theory of Hill and 
Milstein S & §1 and Wang @ are used to discuss branching and the strength of /3-SiC 
under the loading along [001]. 

The present paper is organized as follows. The calculation model is presented in Sec. II, 
where we show the formulation of stress, elastic stiffness coefficients and stability criteria, 
especially the three loading modes with the selection of supercell and the numerical precision 
illustrated at the end of this Section. The benchmark, including equilibrium properties and 
elastic constants, are given in Sec. III. The loading along [111] direction is presented in 
Sec. IV, and biaxial proportional extension is investigated in Sec. V. In Sec. VI, we discuss 
uniaxial extension and tension along [001] direction. Summary and conclusion are given in 
Sec. VII. 

II. FORMULATION 

Consider an unstressed and unstrained configuration, denoted as Xo. It undergoes ho- 
mogeneous deformation under a uniform applied force, and changes from X to X = JX , 
where J is the deformation gradient or the Jacobian matrix. The associated Lagrangian 
strain tensor E is: 

E=i(J r J-I) (1) 
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where T is transpose. The physical strain is: 

e = J - I (2) 

For the present deformation, the internal energy U is a rotational invariant and therefore 
only a function of E. The second Piloa-Kirchhoff stress tensor T is defined as: 

1 dU 

13 ~ V dE~ (3) 
It is related to Cauchy stress, i.e. the true stress by the following equation: 

Tij = det|J| J^ 1 J^au (4) 

where det|J| is the ratio V/V . With the Cauchy stress, the applied force can be obtained 
by multiplying the current transverse area. 

At strained state X, the elastic constants are determined through the following equation: 

C m {X) = —L-^JL^ y =o ) (5) 

where E' is Lagrangian strain around the state X. These elastic constants are rotational 
invariant and symmetric with interchange of indices i «-> j, k <-> / and (ij) <-> (kl), which 
are often expressed in the condensed Voigt notation. 

To analyze the stability, the elastic stiffness coefficient B El is introduced as follows: 

Bijki = Cijki + -{SikO'ji + Sjk<Ju + 5ii<jjk + 5ji<Jik — 25ki<Jij) (6) 

From this definition, we can see that B does not posses (ij) < — > (kl) symmetry generally. 
The system may be unstable when 

det|B| = (7) 

for the first time. 

The following loading modes are considered: 
(i) Uniaxial Extension. 
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i,j = 1,2,3 



(8) 



(ii) Uniaxial Tension. 



i,j = 1,2,3 



(9) 



For a given longitudinal strain, let the transverse lattice contract or dilate to make the total 
energy approach minimum, which corresponds zero stress(traction) on lateral faces. For 
crystal symmetry, the transverse contraction is the same at two perpendicular transverse 
directions, so 



The total energy calculations are carried out with ab initio pseudopotential plane-wave 
program package Fhi96md. @ By means of the mechanism of Hamman @ and Troullier @ , 
the soft first principle pseudopotential S ® is generated, where the local density approx- 
imation(LDA) with the exchange and correlation energy functional developed by Perdew 
and Zunger @ is used. Two supercells are designed in our calculations: one is the 8-atom 
supercell for the equilibrium properties, the loading along [001] and biaxial extension along 
[010] and [001]. The other one is the 6-atom supercell for the loading along [111]; in this case, 
the stacking consequence is Si-C-Si-C-Si-C. There exist two types of {111} plane, between 
Si and C atoms, corresponding to the well known glide and shuffle planes. The glide plane 
cuts three Si-C bonds out of four, and the shuffle plane cuts the remaining Si-C bond. For 
numerical differential feature of stress and elastic constants, the precision must be considered 
carefully. The size of carbon atom is so small that a high cut-off energy is required. Our test 
shows that E cut = 80Ry has also given excellent results.. The k-space mesh is 6 x 6 x 6 for 
the 8-atom supercell and 8x8x4 for the 6-atom supercell in order to keep same precision. 



en = e 22 = -Ae. 
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(10) 



(iii) Biaxial Proportional Extension 



e 2 2 = «e33 /=0 Cij = others 



(11) 
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III. EQUILIBRIUM PROPERTIES 



As the benchmark, we have calculated the lattice constant, elastic constants and moduli 
of /3-SiC of equilibrium. For symmetry of /3-SiC (zincblende structure), there exist three 
independent elastic constants, i.e. Cn, Cu, Cu- The total energy of /5-SiC is calculated 
under the applied hydrostatic, uniaxial deformation, and trigonal strain. Owing to the lattice 
feature of zincblende structure, which includes two fee lattices with a relative displacement 
along [111], the symmetry of center inversion is lost. Four Si-C bonds along directions [111], 
[111], [111] and [111] are not equivalent under the case of trigonal strain. The internal 
atomic position must be fully relaxed, and the internal displacement i3 , which refers to the 
relative displacement of two sublattices beside the displacement from the macroscopic strain, 
will take place. Our results are presented in Table I. From this table, we find, our results 
agree well with the experimental data and the previous first principle and semi-empirical 
calculations. The value of C44 without relaxation, 270 GPa, is already in good agreement 
with the experimental data and better than other theoretical calculations. The relaxed 
value, 254 GPa, is almost the same as the experiment value. Value of the anisotropy A is 
also satisfactory. 

Based on the representation surface il, the moduli along a certain direction can be 
obtained. The Young' modulus along directions [111] and [001] are 554 GPa and 338 GPa, 
respectively. Lambrecht et al. © obtained 603 GPa and 362 GPa. Petrovic et al. @ 
measured Young's modulus of /3-SiC whisker, with an averaged value of 578 GPa with ±10% 
scattering. Applying the equation of the cantilever beam, Wong et. al. HI measured Young's 
modulus of [lll]-oriented SiC nanorod, which are 610GPa and 660 GPa, corresponding to 
the 23.0-nm-diameter and 21.5-nm-diameter SiC NRs. The agreement is good. 

With orientation averages, the moduli of isotropic materials can be obtained. Two 
average methods , namely Reuss averages §1 (E^ and Gr ) and Voigt averages 113 , (Ey and 
G v ) are adopted. According to the theory of Hill 11 , the physical averages, here denoted 
by subscript a, are the intermediate between the Reuss and Voigt averages. With these 
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considerations, Young' modulus and shear modulus of isotropic /3-SiC are given as follows: 

E a ss AASGPa ±2.2% 
G a P3 192GPa ± 2.6% 

Compared with the previous first principle @ and semi-empirical & §1 calculations, our 
results agree better with the experimental values. These results confirm the conclusion of 
Lambrecht: the random orientation hypothesis applies well to the ceramic samples. The 
average Possion ratio u a , 0.17, is close to the experimental value. The small Possion ratio 
of P Si-C, as compared with other materials, e.g. Aluminum 0.347, demonstrates its high 
stiffness . 



IV. LOADING ALONG DIRECTION [111] 

In this loading direction, two loading modes are investigated: uniaxial extension and 
uniaxial tension. At the latter case, the transverse lattices contract in order to approach 
the energy minimum. With regard to the loss of symmetry of center inversion, the four Si-C 
bonds along directions [111], [HI], [111] and [111], are not equivalent under [111] loading. 
The internal relaxation and lateral contraction must be considered. In our calculations, 
by using the lattice constant at room temperature, the internal relaxation is carried out 
after the transverse strains are obtained. The length of Si-C bond is 3.5673 (Bohr) under 
zero loading. As the loading is along [111] direction, the variation of bond length along 
the [HI] direction will be more significant than that of the other three bonds. Fig.l shows 
the energy, stress and force under uniaxial extension and uniaxial tension with or without 
internal displacement. (Without other statement, the strain, force and stress, given in figures, 
are the physical strain, applied force and Cauchy stress.) The corresponding transverse strain 
and internal displacement are given in Fig. 2 and Fig 3. 

At the beginning of loading, as compared with the corresponding Si-C bond length with- 
out the internal relaxation, the difference among the four Si-C bond lengths is small, and 
no marked effect of relaxation and internal displacement through both the shuffle and glide 
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plane is shown. This phenomenon, shown in Fig.l, is obvious for /3-SiC, a kind of high 
stiffness and low Possion ratio covalent material. The strain energy curves of the three 
loading modes, i.e. uniaxial extension and uniaxial tension with or without internal relax- 
ation are almost the same. Despite the fact that the stress and force of uniaxial tension is 
smaller than the uniaxial extension for the relaxation, they are still similar in these loading 
modes. Based on Kleinman's il discussion on Silicon with [111] strain, the internal strain 
tends to keep the bond length along the four unchanged unequivalent [111] directions. In 
our calculations, the internal displacement of atom along [111] direction, which is through 
the shuffle plane, is negative, and that of the remaining three Si-C bonds along the [111], 
[llT] and [Til] directions, through the glide plane, is positive. This means that the effect of 
relaxation always tends to diminish the difference of four Si-C bond lengths. With the in- 
crease of the longitudinal strain, the internal displacement through the shuffle plane becomes 
positive from negative, and that through the glide plane moves from positive to negative 
while approaching zero at the same point ( e z = 0.105). The internal displacements for the 
two cases have the same magnitude but with the opposite sign. The details can be found 
in Fig.3. During the whole uniaxial tension, the magnitude of transverse strain increases 
monotonically. 

With further increase of the longitudinal strain, the strain energy of uniaxial extension 
and of uniaxial tension without internal displacement still have approximately the same 
value. However, both the internal displacements through the shuffle plane and the glide 
plane change their signs ( the symmetry is still hold). The uniaxial tension curve with 
internal relaxation softens quickly, and the shape of tensile curve changes dramatically. An 
energy plateau presents, and material becomes unstable. On the basis of the stress curve, the 
maximum stress of uniaxial tension with internal relaxation, namely the theoretical strength 
a , is obtained and equal to 50.8 GPa. The corresponding critical macroscopic strain and 
internal displacement is 0.144 and 0.082 (Bohr), the Si-C bond length along [111] is 4.163 
(Bohr). With the modified Tersoff potential, Tang and Yip il analyzed the brittle fracture 
of /3-SiC under hydrostatic tension by molecular dynamics. They found that the mode of 
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instability of /3-SiC was the spinodal instability, and the corresponding critical strain and 
pressure were 0.153 and 37.0GPa. Therefore, both the first principle and empirical potential 
calculations gave a similar critical bond-length of /3-SiC along [111]. 

Our result agrees well with the experimental value, given by Wong et al. @ ; 53.4 
GPa, obtained for [lll]-oriented SiC nanorod. This agreement also means that no other 
branching and instability modes exist during the uniaxial loading before it reaches the 
inflexion of energy-strain curve. It is worth while pointing out that the experimental strength 
measured is the bending strength. The tensile and bending strengths are comparable to (3- 
SiC whisker H' El and are also expected to be similar to /3-SiC nanorod. @ Petrovic et 
al. @ measured the tensile strength of /3-SiC whisker and their result is 23.74 GPa, which 
is far smaller than our theoretical calculation and Wong's experiment values for defects. 
Lambrecht et al. @ calculated the tensile strength by Orowan expression with [111] surface 
energy, and the result is 30 GPa. With the similar formula, Op Het Veld and Veldkamp @ 
obtained the theoretical cleavage strength 46.3 GPa, which is close to ours. The detailed 
comparison can be found in Table 2. 

After e z > 0.105, the internal displacement through the shuffle plane becomes positive, 
and that through the glide plane becomes negative. The distance between the atom through 
the shuffle plane along [111] increases, and that through the glide plane decreases, and a crack 
nucleates on the {111} shuffle plane. The internal displacements through the shuffle and 
glide plane at the critical strain are, respectively, 0.082 and -0.082. With further increase 
of the longitudinal strain, the internal displacements of the shuffle and glide planes also 
increase quickly. This positive and negative increase of internal displacements of shuffle and 
glide planes will result in a dramatic cleavage on the {111} shuffle plane and the mixing of 
Si and C atoms through the glide plane. The cleavage on {111} shuffle plane can be partly 
attributed to the lower surface energy than that of the {111} glide plane. 0' El' 10 These 
results agree well with the previous molecular dynamics simulation. © 
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V. BIAXIAL PROPORTIONAL EXTENSION 



To consider only the biaxial proportional extension, this section deals with the extension 
along directions [010] and [001], and not the internal atomic and volume relaxation. The 
strain ratio between [010] and [001] is 0.25, 0.5, 0.75 and 1. The results are shown in Fig.4. 
With the increase of the ratio, the energy, stress and maximum stress will increase at the 
same longitudinal strain accordingly. However, the critical strain is similar for different 
proportional loading modes. 

VI. LOADING ALONG DIRECTION [001] 

In this section, we consider the uniaxial extension and uniaxial tension along direction 
[001]. The reference state is the state with the theoretical lattice constant. Symmetry of 
crystal under this loading mode is tetragonal. Unlike the loading along [111], the deformation 
of the four Si-C bonds in this loading mode is the same and the four bonds are equivalent. 
There will not be any internal displacements and the internal relaxation can be neglected 
during the loading. We have made calculations at several strain with or without internal 
relaxation and found that the value of transverse strain at a specific longitudinal strain is 
the same. Our results are given in Fig. 2, Fig.5 

From Fig. 5a, we can see that both of the strain energy for uniaxial extension and tension 
increase with the increase of the longitudinal strain. The strain energy of uniaxial extension 
is always larger than that of the uniaxial tension, as is expected. However, the energy differ- 
ence between two loading modes is small, same as with the [111] loading. At a larger strain, 
the energy difference becomes even smaller. The applied force and stress of uniaxial tension 
are lower than those of the uniaxial extension for the triaxial stress state at the beginning 
and higher than them at the subsequent loading ( this phenomenon will be explained later). 
Just like Fig. 5a, the difference of applied force and stress between these two loading modes 
is not significant. 
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In order to obtain the ideal strength and analyze the stability under uniaxial tension, 
we calculate the elastic constants and derive the stability criteria based on the stiffness 
coefficients. With the tetragonal symmetry, the number of independent elastic constants is 
reduced to six: C33, C12, C13 = C23, C\\ = C 22 , C44 = C55 and C^; all the other are 
equal to zero. With equation [], [7] and ??, we write the stability criteria as follows: 

(C 33 + <t){Cu + Cia) - 2C 13 (C 13 -a)>0 (12) 

C u - C 12 > (13) 

C 44 + -o- > (14) 

C 66 > (15) 

The first one involves the vanishing of bulk modulus, and is referred as spinodal instability. 
The second instability involves symmetry breaking (bifurcation) with the volume conserva- 
tion; it may be identified as the tetragonal shear breaking and referred as Born instability. 
The third and fourth are two distinct shear deformation instabilities. Six strains are designed 
to calculate the independent elastic constants and given as follows: 



en 


= e 22 = 
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en 


= -e 22 
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In each case, the domain of strain is [0,0.02]. The results are shown in Fig.6 

During uniaxial tension, the variation of mechanical properties on transverse section 
is comparatively small, and the corresponding elastic constants, i.e., Cu, Ci 2 , Cqq, keep 
positive, as shown in Fig. 6a. However, the elastic constants related to longitudinal strain 
change dramatically and even become negative at large strain, e.g. C13 < when 633 > 0.184, 
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C33 < when > 0.352. Because C13 < leads to negative Possion ratio, the transverse 
section will expand with the increase of the longitudinal strain. This phenomenon is also 
shown in Fig. 2, the transverse strain varies en proportionally with 633 when 633 > 0.20. 
Negative Possion ratio had also been investigated by Milstein et al @' 0' & 0. However, 
in their papers, negative Possion ratio only occurs at branching or unstable points, and the 
materials investigated are monatomic metal materials Fe and Ni. In the present calculations, 
it is surprising that /3-SiC, a non-metal and two component crystal, is still stable at this 
negative Possion ratio, shown in Fig.6b. The same result is obtained when calculations 
with internal atomic relaxation are implemented. This phenomenon must be related with 
the bond nature of /3-SiC. The charge transfer §1 and ionic component @ of /3-SiC will 
affect the mechanical response. The detailed analysis of electric structure should be carried 
out and further investigation is necessary. The lattice transverse expansion leads to much 
quicker increase of the force and stress of uniaxial tension than that of uniaxial extension 
and the values of the previous force and stress will be higher than those of uniaxial extension 
at a large strain. 

On the basis of the stability criteria, we have found that the spinodal and Born in- 
stabilities are triggered almost at the same strain 0.37 with the transverse strain -0.0137. 
The corresponding strength, 101.3 GPa, which is almost twice of that for [111]- oriented 
SiC nanorod 53.4 GPa © for its smaller interplannar distance, is obtained. At this criti- 
cal strain, the elongation strain along [111] direction is 0.129, comparable with the critical 
strain 0.144 under uniaxial lading along [111]. Tang and Yip il investigated the instability 
of /3-SiC under hydrostatic tension with stiffness coefficient, and found that failure mode 
of /3-SiC is spinodal instability. This was proved by their molecular dynamics simulations, 
and the nucleation of cracking on the {111} plane and decohesion were revealed. They also 
showed that a shear instability is triggered by the spinodal instability. All of these results 
are same with our first principle calculations. 
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VII. SUMMARY AND CONCLUSION 



On the basis of the DFT total energy calculation and stability theory, we give a detailed 
investigation of mechanical properties of two-atomic constituent materials /3-SiC: energetics, 
elasticity, stress-strain relations, stability and strength under different loading modes and 
directions. The results are satisfactory. 

Owing to the unequivalence of the four Si-C bonds under the uniaxial tension along 
[111], the relaxation must be implemented and internal displacements be calculated. The 
internal displacements along the [111] direction and other three directions namely [Til], 
[111] and [111] have the same magnitude but the opposite sign. At the beginning of loading, 
the effect of relaxation is not significant and tends to diminish the difference of the four 
Si-C bond lengths. However, it becomes important at the subsequent loading and results 
in a crack nucleated on {111} shuffle plane, while the Si atom and C atom through the 
glide plane approach each other. The failure in this loading modes is of cleavage fracture. 
These conclusions are consistent with the previous molecular dynamics simulations. The 
theoretical strength obtained agrees well with the experimental data. 

Under loading along [001], the four Si-C bonds are equivalent for the tetragonal symmetry 
and the relaxation can be neglected. The strain energy, applied force and stress are similar 
despite two distinct loading modes, namely, uniaxial tension and extension. During the 
uniaxial loading along [001], the spinodal and Born instabilities are triggered almost at the 
same strain. Previous molecular dynamics investigation revealed similar facts. Owing to 
the smaller interplannar distance, the corresponding ideal strength, 101 GPa, which is much 
higher than the theoretical strength for the loading along [111] and the experimental data, 
is obtained; however, the Si-C bond length for loading along [001] and [111] at the critical 
strain is close. There exists some stable range when the Poisson ratio is negative, and these 
phenomena are related to the unique bond nature of /5-SiC, a kind of non-typical covalent 
material, which permits charge transfer. A detailed analysis of electronic structure will 
brought up for further investigation. 
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TABLES 

TABLE I. Equilibrium and elastic modulus of /3-SiC. PP-PW, present pseudopotential plane 



wave 


calculations; 


FP-LMTO, Lambrecht (Ref.26); 


CKH, Churcher, Kunc and Heine (Ref. 24); 


Tolpy. 


go (Ref. 22): 


; Tersoff (Ref. 23); 


Exp., expermental values as 


; indicated by footnotes. The 


length unit is bohr 


radius, and modulus is GPa, the anisotropy A= 


2C44/ (Cn 


- C 12 ). 




PP-PW 


FP-LMTO 


CKH 


Tolpygo 


Tersoff 


Exp. 


a Q 


8.166 


8.154 


8.145 




8.164 


8.238 a 


Bo 


225 


223 


224 


211 


220 


225 6 


Cn 


405 


420 




352.3 


420 


390 c 


Cl2 


135 


126 




140 


120 


142 c 


C44 


254(270) 


287 




232 


260 


256° 


A 


1.88 


1.95 




2.20 


1.73 


2.00 c 


Em 


558 


603 




511 


560 


581(±10%) d 610 e 


-EdoO 


338 


362 




272 


367 




Er 


441 


476 




378 


APS) 




Ey 


474 


516 




424 


488 




E a 


458 


496 




401 


475 


448 b 


Gr 


188 


208 




157 


201 




G v 


206 


231 




182 


216 




G a 


197 


219 




169 


208.5 


192 6 


v a 


0.173 


0.146 




0.201 


0.150 


0.168 6 



a Landolt and B6rnstein( |53| ) b Carnahan( p4f) c Obtained from sound velocities ( p5| ) 
d Experimental values from whisker ( |15|) 6 Experemental values from nanorods( [p2]| ) 
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TABLE II. The Young's modulus and strength compared to other theoretical calculations and 
experiments; Here, ts, bs and cs means the tensile strength, the bending strength and the cleavage 
strength. The unit is GPa. 





E 




E" 


a ts 




E c 




E d 




[Ill] 


558 


50.4(0.144) 


603 


30 


46.3 


610 


53.4 


580±10% 


23.74 


[001] 


338 


101(0.37) 


362 















a From FP-LMTO and Orwan expression citelambrecht From Orwan expression |K] 
c Experimental values from nanorods [22] d Experimental values from whiskers |15] 
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FIGURES 

FIG. 1. (a) The calculated strain energy under the uniaxial extension and uniaxial tension 
with or without internal relaxation; (b) Applied force and stress under uniaxial tension with or 
without relaxation. 

FIG. 2. The transverse strain of uniaxial tension along directions [001] and [111]. 

FIG. 3. The internal displacement through shuffle and glide plane under [111] uniaxial tension. 
The club is for the shuffle plane, and the square is for the glide plane. 

FIG. 4. The calculated strain energy(a) and applied stress(b) during biaxial proportional 
extension with difference ratio along [010] (e22) and [001] (633) direction. 

FIG. 5. The calculated strain energy (a), applied force and stress (b) of uniaxial extension and 
uniaxial tension along direction [001]. 

FIG. 6. The calculated elastic constants (a) and stability (b) during uniaxial tension along 
direction [001]. 
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